


VII. Direct and expeditious methods of calculating the excentric 
from the mean anomaly of a planet. By the Reverend Abram 
Robertson, D. D, F.R. S. Savilian Prof essor of Astronomy in 
the University of Oxford, and Radclrffian Observer. Commu- 
nicated by the Right Hon. Sir Joseph Banks, Bart. G. C. B. 

Read February 15, 181 6. 

>ince the publication of Kepler's discoveries in astronomy, 
the attention of men of science has frequently been directed 
to the problem distinguished by his name, and their exertions 
have frequently been employed to overcome the acknowledged 
difficulty of its solution. A statement of the various degrees 
of success, with which these endeavours have been made, is 
foreign to the present design. An account of this kind is now 
also needless, as Dr. Brinkley's examination of such attempts, 
published in the ninth volume of the Transactions of the Royal 
Irish Academy, affords a satisfactory review of most of the 
proceedings on this subject, previous to the year 1802. 

After the following methods had occurred to my consider- 
ation, and I had fully proved their utility by actual application 
to examples, I was anxious to ascertain whether any author 
had anticipated me in the manner in which the investigations 
are conducted. With this view I examined such solutions as 
are referred to in Dr. Brinkley's very able Memoir, all those 
mentioned by Montucla,* of which I could procure a sight, 

# Histoire des Mathematiques, Tom. II. p. 343, &c. I have searched, without 
success, for Lorgna's and Tremblers publications* 
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and some others which had occurred to me in the course of 
my reading, The result of this search is a belief that no one 
before has aimed at a direct solution through the same small 
angle, and, by means of an equation in which this angle and its 
powers are the only unknown quantities, obtained a quickly 
converging series for its value in known terms. The small 
angle being found, with due precision, the eccentric anomaly 
readily becomes known. 

M. Delambre, in his Astronomy, published at Paris in 
1814, in three quarto volumes, calculates the excentric ano- 
maly by a method founded on those of Cassini,* La Caille,^ 
Simpson,^ and Cagnoli.§ This eminent astronomer says of 
it, u Ce procede, le plus directe que je connaises, est au-ssi le 
plus precis ; il n'est qu'appro-ximatif, mais il est toujours 
exact au-dela des dixiemes de seconde pour toutes les planetes 
de notre systeme." That the reader may readily judge how 
far the third method, which I now propose/ deserves atten- 
tion, I have annexed- to my investigation the calculation by it 
of M. Delambre's two examples, and also one relating to the 
comet of 1682 and 1759. 

Each of the following methods of solution is to be considered 
as direct, although it proceeds through the medium of what 
is commonly called Cassini's approximation. This approxi- 
mation, as here used, can only be considered as the first -cer- 
tain step in the computation. No hypothesis is introduced 

* Memoires de P Academic de Pannee 1719. 

f Lecons Elementaires d'Astronomie, Paris, 1761. 

% Essays on several curious and useful subjects, London, 1740.*.. 

§ Trigonometrie, Paris, 1786, 
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into the proceeding, and therefore no correction by trial and 
error is requisite. 




Let ALP be the orbit of a planet, C the centre of the ellipse, 
S that focus in which the sun is placed, and AMP a circle 
described on the greater axis AP as a diameter. Let L be 
the true place of the planet, and AM the corresponding mean 
anomaly. Through L let the straight line ER be drawn, 
perpendicular to AP, and let it meet AP in R and the circle in 
E. Let EC, CM, SM be drawn, and let ST be perpendicular 
to EC, SG to CM and MN to ST. Then MN is parallel to 
ET, and NT is equal to the sine of the arc EM. It is easily 
proved, as in almost every writer on the subject, that ST is 
equal to the arc EM. 

In this Problem it is supposed that AC, CS, AM are given, 
and it is required to find AE the excentric anomaly, for AE 
being found, the true anomaly ASL is easily obtained. 

In each of the three following methods the angles CMS, 
CSM are used, and their difference is found by this pro- 
portion, 

MDCCCXVI. 
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CM + CS t CM— CS : : tan . ± ACM : tan. \ (CSM — CMS), 
Hence the angles become known by their sum and difference* 
As the angle SMN is very small, and consequently the angle 
MSC nearly equal to EGA in the orbits of almost all the 
planets, this way of finding the angle CSM is usually called 
Cassini's approximation to the excentric anomaly ACE. 

First Method. 

Having found the angle CSM, sin. CSM : CM : : sin. SCM : 
SM, which therefore becomes known. Let z equal the angle 
SMN, s equal the series expressing its sine, and c equal the 
series expressing its cosine. Put a equal the sine of CMS,. 
and h equal its cosine. Then, radius being i, ac ™~* hs — s'm» 
CM'N = sin. ECM = sin. (CMS — %). Also, i :SM::5: 
SN = SM x s, and ac — hs + SM x s = TN + SN = EM =-. 
CMS— z, and therefore CMS = % -f ac — bs + SM xs — z + ac 

j^TsMZTb x s. 

Letd=SM — b 9 andthenCMS = z + ac + ds=z + a(i- 

•2.3.4 2-3.4.5.6 l 2.5.4,5.6.7.8 I * V 2.3 * 2.34-5 '' 

. , 7 az* dz 3, . az* . d& 5 

Let ^=CMS— a, and putting A, B, C, &a for the coefficients- 
e — Az- Bz* — Cz 3 + D% 4 + Ez s — F% 6 — G% 7 + &c. 

By reversing this series we find z — — ^ £ + ITTTlblp £ a 

. d 3 ! a* 3 tf 4 , j. 5#J 4 , 5 a 3 4r 

T 6(i + d) 4 * + 2(i+dj** ~~ 24(1 + d) 5 * * I2(i+rf)" 6 ^ ■" 8(i+d)** 

This equation is in parts of the radius, and in order to have 
it in degrees we use this proportion, 1 : 57°, £957795 :: % :5p- 
2957795 % = R ° »> Pitting R° for 57°- $957795* 



z % 
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neiiLuxx <*— ■ ,+^T 2 (i + d) 3 ■ 6(i-h*)* T 2 (i + ( i)s l "^ 

Second Method* 

The substitutions for the sines and cosines of the angles 
CMS, SMN being as in the preceding method, let SG be per- 
pendicular to MC, and then radius being i, i : CS : : sin. 
ACM : CS x sin, ACM= SG. But sin. CMS : SG : : sin. SMN ; 

SN, that is, a : CS x sin. ACM ::s: Sl^l^ s = S N, and 

tlierefore ac — bs -f ?~^h . s ~ XN -j- SN = ST = 

CMS — z. Let ^^r---^^^-— — ~~ b, and then ac <4- d$ — 

Mr " 

CMS — #, and CMS =: z -\- ac -\- ds 9 as in the preceding 
method. 

Third Method, 

Let % = the angle SMN, 5 = the series expressing its sine, 
and c — the series expressing its cosine, as before ; but let a 
now denote the sine of CSM, and b its cosine, and let MN 
meet CS in H. Then ac + bs is equal to the sine of the sum 
of the angles SMH, MSC, that is, ac -\-bs=zz sin. MHA — sin. 
ACE, the excentric anomaly. We have therefore CE : ac + bs 

: : CS : <*££ + *> _ ST = CMS — z. 

Let d~ ~-£™, and*? = ~-~~,and then dc + es^ CMS-— z, 

•» •? dx*" d%/^ dm ^ 

and CMS = % + dc + ^ = % + « h 



2 " 2,3.4 2.3.4.5.6 

+ &c. -f ez — • — + j^j — &c - = « + * + «* — r ~~ zi 

rf «* 4. J2L — &c. Let/ = CMS — d, and then/= As 



2«3>4 2»3"4"*5 

B% 2 *— C% 3 + D^ 4 + E2 5 — &c. putting A, B, C, &c. for the 

coefficients. By reversing this series, or by putting d for a, 
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e for d, and/ for e 9 in the series in the first method, we find 



24(1+ e) s 



+ 12(1 +*) 6 ^ T 8(i+e) ?/ ^ 



I prefer this method to the first or second, and therefore I 
proceed to calculate by it. 

Example L 

Let us suppose with M. Belambre* that the mean ano- 
maly is 135 , and the excentricity of the orbit 0/25, the mean 
distance from the sun being 1. 

HereCM+CS=i.2 5? CM-CS= . 75* and §j^ = ^ 
the log, of which is to be used for any given mean anomaly 



in the orbit. 

Log, tan. 07°, ,30 
Log. tan. 55 .,22,49,844 
CSM = 122 ,.52. 49.844 



~ - Log. 9.7781513 

10.3827757 

10*1609270 



CMS is" found by this 
proportion, 2 062 64^8 ; 1 
: : CMS m seconds : its 
length in parts of the ra-« 
dins. 



CM3 



CS = 0.25 
a 

£^.2099511 
.2115249 



12.. 7m 10.16=43630" 16 

206264,8 

CMS = , 2115249 

Log* 9-39794°o CS 
Log. 9 9241783 b 



Log, 9.3221183 

= CMS 



Lo &* ^ 6 B97$®7 
Log. 5.3144251 

Log, 9.35853616 

L °g- 9*39794oo 
Log, 9.7347108 



.00 1 5738 =/, Log. 7. 19694*95 



£=,1357222 Log. 9.1326508- 
As CSM is obtuse, e is negative. 
1 +£=.8642778 Log. 9.9366534. 



* Vol. II. page 28. 
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The angle SMN is therefore calculated from the series in 
the following manner : 
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First term. 

R° Log. 1.7581226 

/ Log. 7-i9%495 

8-95507 2 i 

1 -f e - - - Log. 9.9366534 

Numb. - .1043322 Log. 9.0184187 
For 2d term .0000231 



Sum 



oM 043553 
SMH. CSM 

ACE = - 



6'..i5."67 

122°„52 ..49, 84 



Second term. 

Log. 1.7581226 

Log. 9.3221183 

- Log. 4.3938990 

15-474*399 

0.3010300 

15; 173 1099 

Log. 9.8099602 
Numb. .0000231 Log." 5-3631497 



RO 

d 

2, Log. 



= - 122..59..5.51 

This differs from M. Delambre's conclusion only in the 
second place of the decimals. 



Example II. 

Supposing, with M. Delambre, that the mean anomaly in 
the same orbit i» g6°, required the excentric anomaly. 

We have as before ^ - - - _ „ Log. 9.7781513 



Log. tan. 
Log. tan. 

CSM = 
CMS = 



•75 
1.25 



48 •• 

33 ■4o , -4i /, .5i 
81 .4,0 4i",5i 
14 "19 •• 1 ^' 49 

CMS 



5155% -49 
206264.8 

.249962 



10.0455626 
- 9-8237 1 39 

Log. 4.7123003 

L °g - 5-gi44 2 5* 
L °g- 9-397 8 75a 
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CS 
a 

d — 



= 0.25 

. 2473S77 

. 2499627 



:CMS 



Log. 9-3979400 
Log. 9-9954030 

Log. 9-393343° 



o 



i 1 + 6 ) 



sd term 



urn 

SMN, 



0025950=/ - Log. 7.4141374 

First term of the series. 

Log. 1.7581226 
Log. 7.4141374 

*Hr« J , l w , w ^ , M ft l ^ W , i n i e % MII , q> t - l , 1 | l , l , 1 l | ,ii b,i 11 - . 

9.1722600 

Log. .0154366 
Log. 9.1568234 



.1434906 
.0000429 

o°. H35335 
CSM 



8'..36".72 

1 ..40 ««4* %5* 



CS = 0.25 



L °g< 9-3973400 

Log. 9.1605669 



= .0361832 Log. 8.5585069 
1+^=1.0361832 Log. .0154366 

Second term of the series. 



R° 


lM:,:,Ml V«/ C»«« • 


1.7581226 


d 
J 


Log. 


9393343° 

4.8282748 


2 - 


Log. 


^•9797404 
.3010300 


(1 +e) 3 


JLOg. 


15,6787104 
.0463098 



Num. .0000429 Log. 5.6324006 



ACE = 81 .49. .18 .23 

This also differs from M. Delambre's conclusion only in 
the second place of decimals. 

Example IIL 
Let us suppose ALP to be the orbit of the comet which 
appeared in 1682, and reappeared in 1759, according to the 
prediction of Dr. Halley ; that CE is equal to 18.07575, that 
that CS is equal to 17.49225, and that the mean anomaly is 
i79°..47'..32 // .i7, it is required to find the excentric anomaly. 

- Log. 8.2149815 



TT prp CM-CS _ .5835 



CM+CS ~~ 

Log. tan. 89°„53 / ..46 // .o8 

83 ..41 ..34 .5 



JL~J\JLlm Idll* 

CSM = 
CMS = 



173 -35 -20 58 
6 ..12 ..11 5 



12.7415531 
a 0.9565346 



= 22331.58 Log. 4.3489194 
206264.8 Log. 5.3144251 
CMS = .1082666 Log. 9.0344943 
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%55 



CS 

a 



1749225 



CE = 18.07575 
d = 



Log. 1,2428457 
Log. 9.0478932 

10.2907389 

Log. 1.2570963 



.1080544 Log. 9.0336426 
.1082666 = CMS 



CS = 17.49225 

b . - 



CE=i8. 07575 



Log. 1.2428457 
Log- 9-9972757 

11.2401214 
Log. 1.2570963 



.0002122 ==/, Log. 6.3267454 
First terra of the series. 



R° 
/ 



(i+0 

Num. .3171788 



Log. 1.758 1226 
Log. 63267454 

8.0848680 
Log. 8.5835638 



R° 

r 
d 



2 



(l+e) 

Num. 



Log. 9.5013042 

Second term of the series. 

Log. 1.7581226 

Log. 2.6534908 

- Log. 9.0336426 

1 3-445' 2 5 6 o 
Log. .3010300 

13.1442260 
Log. 5.7506914 

.0024748 Log. y.s9353^ 

.3171788 

.0000386 

.3196922 sura of the positive 
.0000406 



SMN: 



e = . g6i66y8 Log. 9.9830251 
As CSM is obtuse, e is negative. 

i-f-^.0383322 Log. 8.5835638 
Third terra of the series. 

Log. 1.7581226 
Log. 9.98^0251 
Log. 8.9802362 



R° 

e 

r 

6 



20.7213839 
Log. .7781513 

, , . 4 19-94323 26 

(*+«) - Log. 4.3342552 

Num. — .0000406 Log. 5.6089774 

Fourth term of the series. 

Log. 1.7581226 



R° 

d % 

P - ■ 

(i+o 5 - 

Num. .0000586 



Log. 8.0672853 
Log. 8.9802362 

18.8056440 
Log. .3010300 

18.5046140 
- Log, &.9178190 

Log. 5-5867950 



:o°.3i965i6 = 19'.. 10". 75 

CSM = 173..35..20 .58 



1 73-54 -3 1 S3 = ACE. With the design already 



e — fy — J2L JL J2L ± J^L Sec. But the radius, which is 1, is 

U*/ <y 1 2.3 ' 2.3.4 
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expressed, I adopted in this example the same data with Mr, 
Ivory : see Transactions of the Royal Society of Edinburgh, 
Vol. V. page 236. 

The preceding method bears a nearer resemblance to that 
given by Keill, in his Astronomy, than to any other of which 
I know. Adapting his manner of proceeding to the figure here 
used, he puts y = EM, e = sin. AM,/ cos. AM, and g = CS. 
The series expressing the sine of AE, is therefore equal to 

e t 1 J? 3 » ef 

2 * 2.3 * 2.3. 

to the sine of AE as CS or g is to ST or EM, that is to y, 
Consequently y=ge —gfy ~ ^f- + §£! 4. 1^1 &c. and there- 

fore^ =y-\-gfy + — — ^±L~ J^l!& a 

By reversing this he obtains a series, which, omitting the 
numbers in the coefficients, converges as the powers of -~— > 

CS x sin. AM ^, . , r , . , 

or i + cs x cos. am' 11S ® e g ree °* convergency to the value 

of y in the foregoing examples is as follows. 

In the first example as the powers of -2147372, 
In the second example as the powers of .2553020, 
In the third example as the powers of ,1086700. 
In the third method which has been here investigated, the 

series converges to the value of SMH as the powers of 

/ _ CMS— CS x sin. CSM . r , . , r 

TTe or TTcTT^fcsM • rhls de & ree of convergency m 

the foregoing examples is as follows. 

In the first example as the powers of .0018209, 
In the second example as the powers of .0250438, 
In the third example as the powers of .0055358. 
The third method, by which the three foregoing examples 
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are calculated, appears to me the most simple and precise in 
theory, and the most expeditious in practice of any which I 
have seen. This I say with the greater degree of freedom, 
as I am so well aware of the similarity between its series and 
that of Keili/s, and so perfectly convinced of the advantages 
which it derives from Cassini's approximation, that I consider 
it, with the exception of some deviations, as a combination of 
their methods. 
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